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Abstract. 

The question raised by [Bastin and Martin 2003 J. Phys. B: At. Mol. Opt. 
Phys. 36, 4201] is examined and used to explain in more detail a key point of our 
calculations. They have sought to rebut criticisms raised by us of certain techniques 
used in the calculation of the off-resonance case. It is also explained why this result 
is not a problem for the off-resonance case, but, on the contrary, opens the door to 
a general situation. Their comment is based on a blatant misunderstanding of our 
proposal an d as such is simply wrong. 
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In our paper [2] we have used the Hamiltonian (in the mesa mode case) 

V 2 A 1 
H = 2M + ~2° z + ^ + 2° z) + X ^ z ^ a ~ a + aa+ ^ W 

Let us write equation (1) in the following form (in the mesa mode case f(z) = 1) 
2M 

V = ^-a z + uj(a ] a + ^a z ) + \{a~a) + aa + }. (2) 

It is easy to show that in the 2x2 atomic-photon space the eigenvalues and eigenfunction 
of the interaction Hamiltonian V (V\^) = E^\^)) are given by [3], 



i£ = w(n+i)±^ + A*(n + l), (3) 
|$+) = cos 9 n \n + 1, g) + sin 9 n \n,e), 

|$~) = - sin 6 n \n+ 1, g) + cos 9 n \n,e), (4) 



where 



0„ =tan-' | . AV ^ TT | (5) 

/? +A>+i)-t; 

We write the wave function \^(z,t)) = J2C^(z,t)\^). Then using the total 

n 

Hamiltonian (1) we have 

n ^ ' 

n n 

ri ra 
= E +^<£(*,*)) l^n), (6) 

because of the orthonormality of the wavefunctions |$^) then 



<9t n V 

with no coupling even in the presence of the detuning (equation (13) in BM comment 
[I])- 

It may be worthwhile for the authors to consult some papers that have been 
published previously (see for example Refs. [4-6]) where the detuning has been 
considered and similar results have been obtained. To be more precisely: 
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1- We have used the interaction picture, so that the term (n + does not appear, 
it can be used as a phase only. Bearing in mind the case of mesa mode is being 
treated in our paper i.e. f(z) = 1. 

2- The most serious point is that Bastin and Martin have overlooked the formulae for 

cos2# Tt and sin2# n . (8) 
From equation (10) raised in Bastin and Martin comment, it is easy to write 

t fl AvCT + A 2 (n + l) + f 

tanfc> n = — = -= (9) 

y/f + A2(n + l)-f 

then 

tan2e B = ^I (10) 

2 

Also, it is easy to prove that 

0/3 -a/2 , . ofl __v^tt: 

cos20 n = — ^^^=^^^= and sm20 n = — .(11) 
y/f + A^n + l) ^ + A 2(„ + i) 

Once these formulae inserted in equations (18) and (19) of Bastin and Martin 
comment, we find that, the second terms vanish identically. 

Now let us look more carefully at the general case when we take f(z) no longer a 
constant, i.e. we go beyond the mesa mode case. In this case the orthonormal functions 
m the 2x2 system diagonalize the Hamiltonian V and its elements are diagonal 
in this set of functions with 

\f{z)y/n + l 

tan26» n = s . (12) 

The states |$^) are z-dependent through the trigonometric functions, they satisfy 

dz l n) 1 n) dz' 

l $ n)= ±l $ n)33-| $ n>hjr ■ ( 13 ) 



dz 2 ' nl 1 nl dz 2 1 nl \dz t 
Then \^(z,t)) can be expanded in the form \^(z,t)) = J2^n( z ^)\^n) an d it satisfies 

n 

the Schrodinger equation 

i-jt\*(z,t)) = HMz,t)). (14) 
Hence the coefficients C^(z,t) satisfy the coupled equations 



dz 



.BC- 



dz 
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These equations should replace equations (18) and (19) of the comment of [1]. But once 
f(z) is taken to be constant, then ^ will vanish and we get equations (7) and the results 
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